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Two new approaches to numerical QFT are presented. 



1. Introduction 

"Traditional" Monte Carlo (MC) approaches 
have produced many good results but they do 
have well known limitations. These include re- 
quiring large amounts of computation time as a 
result of the necessity to treat fermions very dif- 
ferently from bosons, as well as problems with 
fermion multiplicity associated with the lattice. 
Moreover, non-positive definite actions, actions 
with rapid oscillations including phase transi- 
tions and symmetry breaking are generally not 
tractable. Finally, highly accurate answers can 
require immense resources. 

For these reasons we have been working on 
two alternative approaches. The first is a tuned 
Monte Carlo method, which we call Mollified 
Monte Carlo (MMC), while the second, involving 
nested approximations to the Schwinger-Dyson 
equations, is cahed the Source Galerkin method 
(SG). These two methods are closely related be- 
cause MMC is initiated by starting from a sta- 
tionary phase point of an action, while the sim- 
plest apphcations of SG tend to favor such points 
in initial iterations. The talks by Ferrante and 
Petrov I^J^ give a few more details. 



portant oscillatory contributions. The resulting 
integrals are evaluated around stationary points 
using appropriate importance function weighting 
and MC integration. In the cases we have exam- 
ined this approach allows accurate evaluation of 
normally intractable integrands, including com- 
plex integrands. In particular, theories can be 
directly evaluated in symmetry breaking or other 
phases with comparative ease. MMC still suf- 
fers from the traditional problems associated with 
fermions (although with the potential for more 
rapid convergence of integrals) but does consid- 
erably extend the class of actions and regions of 
evaluation open to the usual approaches. 

3. Source Galerkin Methods 

Our Source Galerkin method, while requiring 
less familiar numerical techniques, is powerful be- 
cause it is defined on the continuum and treats 
fermions, except for anti-commutativity, in the 
same way as bosons. The basic idea is simple: 
We assume that the QFT action is written with 
sources Ji (x) for every field (x) so that the am- 
plitude Z = (+0|0— ),/. satisfies the differential 
equation. 



2. Mollified Monte Carlo 

Unlike normal MC, the integrands for problems 
approached using MMC are very carefully con- 
trolled before any actual heavy duty calculation 
takes place. Regions of high oscillation are iden- 
tified and smoothed (mollified) by using simple 
exponential functions that do not change the val- 
ues of any of the integrals but do suppress unim- 
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essentially the 
field equations 



A familiar but far from straightforward exam- 
ple is (j)^ scalar QFT, 

{d^ + m^) 4){x) + g 4)^{x) = J{x) , 
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which becomes (in Euclidean space), 
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Z[J] = 



This is a very non-trivial infinite set of coupled 
differential equations. To get some idea of its 
complexity, examine the one point (0-dim) case. 



d 



2d I 



Z\J\ = , 



Even this simple equation has three independent 
solutions and requires the input of parameters for 
a full solution |Q. These three solutions can be 
characterized by their degree of singularity: 

• Regular at g — > 0: consistent with pertur- 
bation theory; 

• ^ <?~^^^ at g — > 0: "symmetry breaking"; 

• ~ exp(^^/4(7) at g — > 0: "instanton". 

In the path integral language, these solutions cor- 
respond to the three saddle points of the path 
integral in the complex 0-plane. 
The Schwinger condition, 

J {dx] { -dr^ Z 
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dg 



is used to stabilize the phase choice |^ using SG in 
iterative approximations as is discussed in what 
follows. 

A direct approach to solving these is to assume 
a truncated functional expansion in the sources 
and introduce some way to minimize the error. 
It is the presence of a method for error control 
that differentiates our approach from any older 
Schwinger-Dyson approximations. Formally, the 
correct answer is, 

Z = cxpi^Ao+j Ai{x) Ji{x)+j A,j{x,y) Ji{x)Jj{y)+- ■ ■ | 

The problem is to introduce a procedure to iter- 
atively calculate the AiS. In order to do this, we 
put information we already have about the field 
theory (spectral information, relativity, Lehman- 
Kallen representations) and make decompositions 
by assuming 

Z = expja + Jn J(a::) + J J{x) G{x - y) J(y)| 

+ {all possible higher order Green's function terms 
all constructed from terms involving only functions 
G{x — y) but with possibly different parameterization.} 



The above representation is highly symboHc. 
To make calculations tractable we have spanned 
the solution space involving products of sources 
by multiplying these, in turn, with all possible 
products of free field propagators with arbitrary 
weights and masses. For example, the quartic 
source contribution, JJG4JJ is represented by all 
possible graphs involving only the "free" G{x — y) 
with four external lines. The G{x—y) can have ar- 
bitrary numerical weights and mass. The result is 
a very general, very complicated expression con- 
sistent with Lorentz and translation invariance. 
The propagator G{x — y) is. 



G{x - y) = jdK^d^k- 
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This satisfies space-time restrictions, and the 
a{K'^) are weight functions. When inserted into 
the field equations the above expansion yields 
constraints on a{K'^). In general, this is too com- 
plicated to directly solve, since this would be tan- 
tamount to an exact solution of a field theory. Fi- 
nally, there are cross-terms in any equation and 
the equations are non-linear! 
We simplify as follows: 

• Truncate the expansion in J; 

• Limit the number of masses in each propa- 
gator; 

• Limit the number of graphs considered for 
each J. 

We do this in an organized, systematic, way so 
that after a first guess more terms can be included 
so as to iterate answers. 

Roughly this works as follows: Approximating 
yields Z*[J] and, of course. 
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[J]^0. 



The idea of Source Galerkin is to require that, 

J dJi ■ ■ ■ dJn Fi{Ji, . . . , J„) FyTJI Zapprox[>^] = 



\6J 

so that Zapproxi'^] Satisfies the field equations "on 
the average". An appropriate number of Fi are 
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picked so that all the undetermined weights in 
Z*[J] are determined. 

In general, these equations are non-linear, and 
must be solved in a very careful and system- 
atic manner. Theorems for simpler Galerkin ap- 
proaches promise convergence and we conjecture 
that the same is true for QFT. We have applied 
this approach to many models (with amazing ac- 
curacy when a check is available) but the stability 
and convergence must be confirmed. Higher iter- 
ations are simple in principle, but introduce com- 
putational difficulties. To understand these we 
have re-examined perturbation theory using vari- 
ants of our numerical techniques [4] and also re- 
analyzed trivial 0-dimensional models for higher 
order Galerkin expansion [2]. 

4. Perturbation Theory 

The usual graphical rules apply (not our ex- 
tended "exact" rules). Application of our ideas 
have lead to a new way to numerically calculate 
graphs 1^. The major ingredients of our new pro- 
cedure are the Sinc-function expansion: 



Sk{h,x] 



_ sin(7r [x — kh)/h) 
■K (x — kh)/h 
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and the cutoff propagator 



Ga(x) 
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Using the Sine expansion gives [4] 
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Typically, this can be approximated to very 
high accuracy (1 part in 10^^) with fewer than 100 
terms in the sum and similar statements hold for 
fermion propagators. Using these propagators, 
graphs can be reduced to multi-dimensional sums 
which can then be quickly (relative to the anal- 
ogous Monte Carlo integral) calculated to high 
precision. We believe this may be a powerful 



way to check high order magnetic moment cal- 
culations, which are currently of great interest. 
Technical issues with accuracy using an "auto" 
renormalization based on approaches used in lat- 
tice gauge theory scheme have slowed us down, 
but we hope to release some more computations 
emphasizing the power of this method in the near 
future. Moreover, the ease and speed of these per- 
turbative calculations gives us hope that we can 
iterate Source Galerkin calculations to moderate 
order with relatively small amounts of computer 
time. 
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